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Neutron stars are believed to contain (neutron and proton) superfluids. I will 
give a summary of a macroscopic description of the interior of neutron stars, in a 
formulation which is general relativistic. I will also present recent results on the 
oscillations of neutron stars, with superfluidity explicitly taken into account, which 
leads in particular to the existence of a new class of modes. 

1 Superfluid relativistic hydrodynamics 

Neutrons stars are believed to contain in their interior superfluid neutrons and 
superconducting protons. This is suggested by nuclear physics calculations and 
by the very long relaxation time scales after 'glitches' (sudden increases of the 
neutron star angular velocity). Over the last years, we have been developing a 
formalism that can describe the relativistic hydradiiaamics of mixtures of superfluid 
or superconductors and possibly ordinary fluids Id'cHj. 

The hydrodynamics of several species labelled by the index X, with respective 
partiele currents and respective electric charge per partiele , can be described 
by the Lagrangian 

C = Am + J2e-nPAp + ^Fp^F-'P , (1) 

where Ap is the electromagnetic gauge form and Fpa the corresponding electro- 
magnetic tensor, and where Am is the 'hydrodynamical' part of the Lagrangian, 
depending only on the partiele currents n^, the variations of which define the mo- 
mentum covectors: SAm — fJ-^p^nP. This variational principle, corresponding to 
perfectly conducting fluids, leads to matter conservation equations, one for each 
species, Vp rif — O (one can generalize to allow for chemical reactions between var- 
ious species □) , and to Euler-type equations of motion, which can be written in the 
very compact form 

K^IpK]^^^ ^itli TT- = /i-p + e-Ap (2) 

In order to deal with superfluids or superconductors, one must impose the fur- 
ther condition that the momentum covectors ir^p are gradients. This is in fact true 
on small scales only, because on larger scales the superfluid or the superconductor 
will be threaded in general by arrays of vortices, due to a global angular momen- 
tum (for superfluids) or an external magnetic field (for superconductors). The 
macroscopic description, which takes into account the average effect of the vortices, 
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can be done by considering the (generalized) vorticity tensors w^^ = 2V[p7r^j, now 
non-zero since there are vorticies, as fundamental quantities i. 

2 Oscillations of superfluid neutron stars 

In a neutron star, the existence of a superfluid component, weakly connected to 
the normal component, leads to a richer spectrum of oscillations. In general rela- 
tivity, one must take into account not only the matter perturbations but also the 
perturbations of the metric 5^,^. In a specific gauge, the even-parity modes can be 
written (the study of the m = O modes is enough because of the degeneracy in m) 
in the form 

Sgrr = -eV'ff2e*"*Pi(Ö), ögee = Sg^^/ sin" 9 = -r'+^if e*"*Pz(ö), (3) 

where Pi stands for the Legendre polynomial of order L In a two-component model, 
the matter perturbations can be described by the two matter displacements ^„ and 

C - r'-ie-V2w„e^-*Pi(ö), C = -r'-^Vry^'dePi{e), 

Cp = r'-^e-^^^Wpe'^'Pii9), ^ = ^r'-^Vpé'^'dePm. (4) 

Inserting the above expressions for the perturbations into the perturbed Einstein's 
equations as well as the perturbed Euler equations, one ends up with a system of 
linear equations, consisting of two contraints (one is H2 = Hq, the other expresses 
Hq in terms of the other perturbations) and of first order differential equations of 
the form 

^ = Qi^^Y^ (5) 

where Y is 6-dimensional column vector containing Hi, K, Wn, Vn, Wp and Vp, 
and Qi_cü is a 6 x 6 matrix with r-dependent coefiicients that depend only on the 
background configuration, as well as on l and co. Considering the boundary con- 
ditions at the center and at the surface of the star, this system can be solved, up 
to a global amplitude, for any value of uj. The physically relevant modes, however, 
also called quasi-normal modes, correspond the specific values of lu for which the 
metric outside the star represents only outgoing gravitational waves. 

A numerical investigation for a very crude model of two independent poly- 
tropes (with different adiabatic indices) has shown that the two-component star 
will exhibit new modes, superfluid modes, which are specific of the existence of two 
components since the two fluids are counter-moving for these modes □. 
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